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ABSTRACT 


Previous  works  by  others  have  presented  evidence  of  a  satellite 
resonance  phenomenon  with  geopotential  tesseral  terms  of  order  m  =  13  and 
m  14.  This  report  presents  an  analysis  of  the  data  of  five  satellites  and 
provides  evidence  for  this  resonance  phenomenon  with  terms  of  order  m  =  16, 
m  17,  and  m  32.  By  using  two  atmospheric  density  models  in  the  analy¬ 
sis,  it  was  concluded  that  even  in  a  high-drag  environment  one  observes  an 
underlying  satellite  resonance  phenomenon  with  geopotential  forces  arising 
from  high-degree  and  high-order  tesseral  terms  of  the  geopotential  expan¬ 
sion.  An  accurate  determination  of  the  descriptive  parameters,  difficult  with 
the  present  low -altitude  satellite  data  and  atmospheric  density  models  avail¬ 
able,  must  therefore  be  the  result  of  a  joint  geopotential-atmosphere  study, 
or  await  the  direct  use  of  drag  accelerometer  data. 
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the  product  of  the  earth's  mass  and  the  constant  of  gravitation  G 
the  earth 's  radius 

Legendre  polynomials  of  the  first  kind 
the  associated  Legendre  polynomials  of  the  first  kind 
the  geocentric  radius 
tlu1  latitude 

Greenwich  oast  longitude 


Results  presented  here  are  derived  primarily  from  an  analysis  of  the 
tracking  data  of  satellites  1962  BE  1,  1963  55A,  and  1965  79A,  the  nominal 
orbital  characteristics  of  which  are  presented  in  Table  I.  The  raw  data  are 
range,  azimuth,  and  elevation  obtained  by  the  tracking  network  of  the  Air 
Force  Satellite  Control  Facility.  ; 

Table  I  also  contains  nominal  orbital  data  for  the  considerably  higher- 
perigee  satellites  1963  35A  and  1963  42A.  All  the  results  presented  in  this 
report  for  satellites  1963  35A  and  1963  42A  are  the  work,  respectively,  of 
Mr.  Lem  Wong  and  Mr.  Richard  J.  Farrar,  both  of  Aerospace  Corporation, 

Table  I.  Nominal  Satellite  Orbital  Characteristics 


Satellite 

Perigee 
Height 
(n  mi) 

Eccentricity 

Inclination 

(deg) 

Nodal 

Period 

(min) 

1962  BE  1 

125 

0.  017 

82.  0 

91  0 

1963  55A 

105 

0.  014 

64.  9 

90.  0 

1965  79A 

115 

0.  009 

75.  1 

89.  7 

1963  35A 

165 

0.  002 

81.  9 

90.  8 

1963  42A 

1  55 

0.  005 

89.  9 

90.  9 

SECTION  II 


PROCEDURES 


The  first  step  was  to  produce,  for  each  pass  of  raw  data  reported,  a 
smoothed  reduced  data  set  that  represents  the  rrore  numerous  raw  data  of 
the  pass.  This  was  done  by  making  a  local  least  -  square  s  fit  for  each  pass 
of  data  and  computing  three  sets  of  range,  azimuth,  and  elevation  data  from 
the  resulting  best-fit  6-vector  of  position  and  velocity  at  the  local  epoch. 

The  remaining  and  principal  portion  of  the  study,  now  to  be  described,  was 
accomplished  by  using  TRACE,  the  Aerospace  Corporation  Orbit  Determina¬ 
tion  Computer  Program  basically  described  in  Reference  3. 

Procedures  using  TRACE  are  primarily  of  two  types  and  can  be  out¬ 
lined  by  the  four  following  steps: 

a.  Through  long-arc  fitting  to  the  compacted  data,  a  best-fit 
7-vector  of  position  and  velocity  at  epoch,  and  the  drag 
parameter  CpA/W  are  obtained  (where  Cp  is  the  drag 
coefficient,  A  the  frontal  area,  and  W  the  weight  of  the 
satellite).  For  the  near-earth  satellites  investigated, 
geopotential  forces  and  atmospheric  drag  have  been 
considered  sufficient  representations  of  the  forces  acting 
upon  these  satellites.  The  basic  geopotential  used  through¬ 
out  this  study  has  been  a  6th-degree  and  order  model  of 
Guier,  which  is  unavailable  in  the  open  literature.  Some 
experimentation  with  atmospheric  density  models  was  made 
by  using  the  Lockheed- Jacchia  (Reference  4)  and  Jacchia 
1964  (Reference  5)  models. 

b.  In  order  to  determine  how  well  the  best-fit  7-vector,  and 
the  force  model  used  in  its  determination,  represent  the 
motion  of  the  satellites,  the  next  step  is  to  hold  these 
quantities  fixed  and  to  solve  for  noncor related  time  biases 

for  each  pass  of  data.  In  the  absence  of  any  true  clock  errors, 
these  time  biases  may  be  interpreted  as  errors  in  prediction 
of  the  time  of  arrival  of  the  satellite  at  a  station  due  to  vector 
and  force  model  errors,  station  location  errors,  etc.  One  may 
obtain  in-track  position  error  estimates  by  multiplying  the  time 
biases  by  minus  inertial  space  velocity  (i.e.  ,  speed)  of  the 
safellite.  For  satellites  considered  here,  a  nominal  value  for 


-3- 


this  speed  is  25  feet  per  millisecond.  The  minus  (speed) 
results  from  the  fact  that  a  negative  time  bias  signifies  a 
predicted  arrival  at  the  station,  whi~h  is  earlier  than  what 
is  actually  observed,  and  this,  in  turn,  signifies  a  positive 
in-track  position  prediction  error. 

c.  If  examination  of  the  time  bias  patterns  obtained  in  Step  b 
reveals  any  gross  deterministic  force  model  errors,  one 
returns  to  Step  a  and  augments  the  state  vector  to  include 
parameters  that  describe  the  candidate  force  model  error 
selected. 

d.  In  order  to  detern  <ne  how  well  the  candidate  force  model 
error  removes  the  1.1-track  error  presented  in  the  time 
bias  pattern  of  Step  b,  one  returns  to  Step  b  with  the  best-fit 
augmented  state  vector  of  Step  c. 

In  Figures  1  through  7,  are  presented  the  time  bias  patterns  resulting 
from  the  best-fit  7-vector  for  satellites  1962  BE  1,  1963  55A.  1965  79A,  and 
1963  42A;  those  for  1963  35A  were  not  available  in  the  present  fojmat.  The 
force  model  used  is  indicated  in  each  figure.  In  all  cases,  we  observe 
periodic  variations  in  the  time  bias  patterns  similar  in  nature  to  the  variations 
in  the  along-track  residual  patterns  of  References  1  and  2. 

The  first  7-vector  time  bias  patterns  obtained  were  those  employing  the 
Lockheed- Jacchia  atmospheric  density  model.  However,  since  most  of  the 
satellites  studied  are  in  a  relatively  high-drag  environment,  we  turned  to  the 
Jacchia  1964  model  (the  most  sophisticated  model  locally  available)  to  see  if 
these  variations  might  not  be  removed  by  a  supposedly  improved  density 
model.  The  results  of  this  atmospheric  density  model  experiment  may  be 
seen  by  examining  Figures  1  through  6. 

The  most  important  fact  revealed  is  that  even  though  the  amplitude  of 
the  time  bias  patterns  sometimes  shows  a  strong  dependence  on  the  atmos¬ 
pheric  density  model  used,  the  period  of  these  variations  is  relatively 
insensitive  to  the  choice  of  the  density  model. 

It  may  still  be  argued  that  these  variations  in  the  time  bias  patterns 
are  due  to  some  remaining  error  in  modeling  the  atmospheric  density 
variations  experienced  by  the  satellites  (o-r,  more  generally,  to  errors  in 
modeling  the  drag  acceleration);  the  principal  effort  of  this  study,  however. 
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Figure  1.  T  mie  Bias  Pattern  without  Resonant  Terms:  Satellite  1962  BE 
( Lockheed -Jacchia  Atmospheric  Density  Model) 
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Figure  2.  Time  Bias  Pattern  without  Resonant  Terms:  Satellite  1962  BE  1 
(Jacchia  1964  Atmospheric  Density  Model) 
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without  Resonant  Terms:  Satellite  1963  55A 
Atmospheric  Density  Model) 
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Figure  5.  Time  Bias  Pattern  without  Resonant  Terms:  Satellite  1965  79A 
(Lockheed-Jacchia  Atmospheric  Density  Model) 
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Figure  6.  Time  Bias  Pattern  without  Resonant  Terms:  Satellite  1965  79A 
(Jacchia  1964  Atmospheric  Density  Model) 
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Figure  7.  Time  Bias  Pattern  without  Resonant  Terms:  Satellite  1963  4?A 
( Lockheed -Jacchia  Atmospheric  Density  Model) 


is  to  show  the  extent  to  which  these  observed  variations  in  the  time  bias 

pattern  may  be  attributed  to  resonant  tesseral  terms  omitted  in  the 

geopotential  expansion.  Having  an  observed  value  for  the  period  of  the 

resonance  phenomenon  from  the  time  bias  patterns,  we  next  consider  the 

theoretical  prediction  for  likely  geopotential  terms. 

From  Reference  2,  one  may  obtain  the  theoretically  expected  period 

of  the  resonance  phenomenon  PD  as 

K 


p  =  Ztt/U 
R  [j  +  mp  | 


(2) 


where  m  is  the  order  of  the  geopotential  term  j  =  1,  2,  3,  •  •  ■  ,  and 
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for  the  satellites  discussed  here. 

In  Eq.  (3) 

oo  =  the  earth's  sidereal  rotation  rate 
e 

n  =  Keplerian  mean  motion  of  the  satellite 
oo  =  argument  of  perigee 

ft  =  inertial  longitude  (right  ascension)  of  the  ascending  node 

The  mean  argument  of  latitude  U  is  equal  to  tne  sum  of  oo  and  the 
mean  anomaly  M. 

As  a  very  rough  approximation,  one  may  take 
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and  generate  curves  such  as  those  presented  in  Figure  8.  One  may  then  use 
Figure  8  to  obtain  a  rather  good  qualitative  estimate  for  nearness  toresonanc 
simply  from  the  Keplerian  per’od  P^. 

A  more  accurate  determination  of  can  be  obtained  by  using  the 
first-order  oblateness  expressions 


and 


(6) 


(7) 


(8) 


in  Eqs.  (2)  and  (3),  where 

=  the  coefficient  of  the  second  zonal  harmonic 
a  =  the  major  semi-axis 
e  =  the  eccentricity 

i  =  the  inclination  of  the  satellite  orbit. 

Earlier  a^  was  defined.  From  the  foregoing  equations,  one  may 
construct  tables  of  for  a  given  satellite  as  a  function  of  j  and  m. 

Table  II  presents  such  data  for  satellites  1962  BP!  1,  1963  56A,  and  1965  79A. 
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Table  II.  Theoretical  Values  of  the  Period  of  Resonance  PR  Phenomenon 


Satellite 

j 

m 

PR  (days) 

1962  BE  1 

1 

1  5 

1 . 4 

1 

16 

3.  3 

1 

17 

0.  8 

2 

31 

2.  4 

2 

32 

1.  7 

2 

33 

0.  6 

1963  55A 

1 

1  5 

1.  3 

1 

16 

4.  5 

1 

17 

0.  8 

2 

31 

1.8 

2 

32 

2.  2 

2 

33 

0.  7 

1965  79A 

1 

15 

1.1 

1 

16 

8.  6 

1 

17 

0.  9 

2 

31 

1.  3 

2 

32 

4.  3 

2 

33 

0.  8 

As  an  example  of  the  use  of  Table  II,  a  comparison  of  the  3.  1-day 
period  in  the  time  bias  patterns  for  satellite  1962  BE  1  (see  Figures  1  and  2) 
with  the  entries  ir.  Table  II  strongly  suggests  that  the  missing  resonant  term 
would  have  order  m  =  16  (and  j  =  1).  By  returning  to  the  theory  presented  in 
Reference  2,  it  can  be  seen  that  odd  j  signifies  odd  degree  n,  with  n  =  17 
being  the  first  odd  degree  for  m  =  16.  Res.  Its  for  the  tesseral  term  (n,m)  = 
(17,  16)  on  satellite  1962  BE  1  will  be  presented  in  the  next  section,  along 
with  those  for  the  remainder  of  the  flights  in  Table  I.  Table  III  presents  a 
comparison  of  the  observed  and  predicted  resonance  periods  for  satellites 
1962  BE  1,  1963  55A,  and  1965  79A,  and  the  associated  resonant  term 
predicted.  Similar  small-percent  discrepancies  between  the  predicted  and 
observed  resonance  periods  noted  in  Table  III  were  also  reported  in 
Reference  1. 
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Table  III.  Comparison  of  the  Observed  and  Predicted  Resonance 
Periods  and  Associated  Resonance  Term 


SECTION  III 


RESULTS  AND  DISCUSSION 


Table  IV  presents  a  summary  of  the  results  on  resonance  terms  derived 
from  the  satellites  listed  in  Table  I.  Figures  9  through  16  show  the  time  bias 
patterns  for  satellites  1962  BE  1,  1963  55A,  1965  79A,  and  1963  42A  that 
prevailed  when  the  derived  values  in  Table  IV  were  used  in  the  force  model. 

To  determine  how  well  the  derived  terms  do  remove  the  in-track  error 
observed  earlier,  Figures  9  through  16  should  be  compared  with  the  corre¬ 
sponding  7-vector  solutions  presented  in  Figures  1  through  7. 

In  ail  combinations  of  the  satellites  and  atmospheric  density  models 
tested,  one  observes  a  nearly  complete  removal  of  all  systematic  variations 
in  the  time  bias  patterns  of  Figures  1  through  7  by  the  derived  resonant  term 
predicted  for  the  individual  satellites.  Further,  Table  IV  shows  that  the 
derived  coefficients  of  the  same  degree  and  order  are  generally  in  close 
agreement  with  each  other  and  with  rule-of-thumb  values  suggested  by  Kaula. 

Regarding  this  last  statement,  in  a  private  communication  Kaula  has 

suggested  that  an  examination  of  autocovariance  analysis  values  indicates 

that  E{Jnm}  =  0.  07  X  10  for  a  degree  n  in  the  range  13  to  19.  Unnormaliz- 

mg  this  value  for  starting  values  for  (n,  m)  =  (17,  16)  and  (n,  m)  =  (17,  17) 

-  25  -  26 

yields,  respectively,  2.  0  X  10  and  3.  4  X  10  ,  which  are  in  generally 

close  agreement  with  the  values  derived.  This  is  particularly  true  of  the 
highe r -altitude  satellites  1963  35A  and  1963  42A,  on  which  the  effect  of  the 
atmosphere  is  less  and  there  is  less  variation  in  results  when  changing 
atmospheric  density  models. 

The  discordant  value  obtained  for  (n,  m)  =  (17,  16)  from  satellite  1962 
BE  1  when  using  the  Jacchia  1964  atmospheric  density  model  has  subsequently 
been  shown  to  be  due  to  the  manner  in  which  this  model  incorporates  density 


variations  correlated  with  the  earth's  magnetic  activity.  3y  setting  the 

planetary  magnetic  amplitude  at  a  constant  value  in  the  Jacchia  1964  model, 

-  25 

one  obtains  ^  =  3.  29  X  10  and  ^  =  -13.  58  degrees  from  the  data 
of  satellite  1962  BE  1. 
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BASIC  GEOPOTENTIAL  MODEL.  GIIIER  6th  DEGREE  AND  ORDER  (UNPUBLISHED) 
RESONANT  TERM  (SOLVED  FOR):  JI7,  IS  0.42940699E-24  LI7,  16  -13  4  51558  DEG 


(Lockheed -Jacchia  Model) 


BASIC  GEOPOTENTIAL  MODEL  SIXER  6  th  DECREE  AND  ORDER  (UNPUBLISHED) 
RESONANT  TERM  (SOLVED  FOR):  JI7,  16  0. I5426335E-23  LI7,  16  -14.165219  DE6. 


BASIC  GEOPOTENTIAL  MODEL  GUIER  6th  DEGREE  AND  ORDER  (UNPUBLISHED) 
RESONANT  TERM  (SOLVED  FOR):  JIT,  16  0.49927375E-24  LI  7,16  -13.506150  OEG 
RESONANT  TERM  (SOLVED  FOR):  J32, 32  0.1 8024650E- 50  L32.32  3.8267676  DEG 


Figure  11.  Time  Bias  Pattern  with  Resonant  Terms:  Satellite 
(.Lockheed -Jacchia  Model) 


Time  Bias  Pattern  with  Resonant  Terms:  Satellite  1963  55A 
(Lockheed-Jacchia  Atmospheric  Density  Model) 


BASIC  GEOPOTENTIAL  MODEL  GOIER  6th  DEGREE  AMD  ORDER  (UNPUBLISHED I 
RESONANT  TERM  (SOLVED  FORI:  JI7.I6  0.3244J987E -24  LI7, IS  -  7 1288285  DEG 


Figure  13.  Time  Bias  Pattern  with  Resonant  Terms:  Satellite  1963  55A 
(.lacchia  1964  Atmospheric  Density  Model) 


BASIC  GEOPOTENTIAL  MODEL  :  tUIER  61*  DEGREE  ANO  ORDER  (UNPUBLISHED! 
RESONANT  TERM  (SOLVED  FOR) :  J)2,  32  0.42B82I65E-50  L32  ,  32  +0.079735 DEG 


Figure  14.  Time  Bias  Pattern  with  Resonant  Terms:  Satellite  1965  79A 
( Lockheed -Jacchia  Atmospheric  Density  Model) 


BASIC  SEOBOTERTIAl  MODE  l  :  ei':tR  cm  OCMCE  AMO  OR  DEB  (UNRVW.ISMEO) 

RESONANT  TERM  (SOLVED  FOR) :  J32,  32  0.S2S2S624E  -  SO  132,32  -I.32S2SI  DCS 
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Figure  15.  Time  Bias  Pattern  with  Resonant  Terms:  Satellite  1965  79A. 
(Jacchia  1964  Atmospheric  Density  M^del) 


BASIC  GEOPOTENTIAL  #0 DEL  GUIER  6t»  DEGREE  AND  ORDER  (UNPUBLISHEDI 
RESONANT  TERM  (SOLVED  FOR)  .'17,16  0  294542J9E-24  LIT,  16  -I0.5M57I  Dt( 
J 1 7, 17  0  J65842I6E-25  U7.I7  -HJ4SIS26  DEI 
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igure  16.  Time  Bias  Pattern  with  Resonant  Terms:  Satellite  1963  42A 

(Lockheed-Jacchia  Atmospheric  Density  Model) 


SECTION  IV 


CONCLUSION 


Atmospheric  density  variations  improperly  modeled  may  contribute  to 
the  time  bias  patterns  of  Figures  1  through  7,  but  neither  the  Lockheed- 
Jacchia  nor  the  Jacchia  1964  models  represent  perfectly  the  density  variations 
experienced  by  these  satellites.  By  recognizing  these  two  facts  simultaneously, 
it  is  concluded  from  the  foregoing  results  that  even  in  a  high-drag  environment 
one  still  observes  an  underlying  satellite  resonance  phenomenon  with  geo¬ 
potential-forces  arising  from  high-degree  and  order  tesseral  terms  of  the 
geopotential  expansion. 

It  is  recognized,  however,  that  an  accurate  determination  of  the  descrip¬ 
tive  parameters  would  be  extremely  difficult  from  the  present  low-altitude 
satellite  data,  especially  with  the  atmospheric  density  models  presently 
available.  Their  determination  must  result  from  a  joint  geopotential- 
atmosphere  study,  which  will  refine  both  members  of  the  force  model. 
Otherwise  the  more  exotic  techniques  of  drag-free  satellites  (i.  e,  ,  no  drag 
force  modeling  required)  or  the  direct  use  of  drag  accelerometer  data  must 
be  awaited. 
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